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I review the physical and mathematical foundations for the theoretical description of 
the hadron wavefunction at small a; as a Colour Glass Condensate. In this context, I dis- 
cuss the phenomenon of gluon saturation and some of its remarkable consequences: a new 
^ ! "geometric scaling" for F2, which has been recently identified at HERA, and the unita- 
\Q \ rization of the hadronic cross-sections at high energy. I show that by combining saturation 
^ ' and confinement one obtains cross-sections which saturate the Froissart bound. 
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1. INTRODUCTION 

Understanding the high-energy behaviour of hadronic cross-sections is a most intriguing 
and fascinating problem, which is intimately related to the physics of high parton densities. 
, Linear evolution equations in QCD, like DGLAP and BFKL, predict a rapid growth of the 
Qj [ parton densities, which violates unitarity constraints. This is the "small-x problem" Q of 
'-^ perturbative QCD. But the applicability of perturbation theory to high-energy problems 
is by no means obvious. Indeed, even when s ^ 00 (with s the total energy squared in the 
^ ■ center of mass frame), there is a priori no guarantee that all the transferred momenta are 
large. The "infrared diffusion" problem of the BFKL equation [ Q (see Sect. 2.1 below) 
provides an explicit counterexample in this respect. 

But this does not prove the inadequacy of perturbation theory in general. In fact, the 
strong rise of the parton densities carries in itself the potential for its own solution [ |^, ^, 
^ 1^ : The interactions between the partons from different parton cascades, omitted in the 
linear evolution equations, will become more and more important with increasing energy, 
and may eventually tame the growth of the parton densities. This parton saturation 
phenomenon is expected to introduce a characteristic momentum scale, the saturation 
momentum Qs{x, A), which is a measure of the density of the saturated gluons, and grows 
rapidly with 1/x and A (the atomic number). Thus, for sufficiently high energies and/or 
very large nuclei, gluon saturation provides an intrinsic hard momentum scale which limits 
the "infrared diffusion" and justifies the use of weak coupling methods: as{Qs) ^ 1- 

Note, however, that although the coupling is weak, the saturation regime remains non- 
perturbative in that the non-linear effects associated with the high parton densities cannot 
be expanded out in some perturbative scheme, but must be included exactly. This is a 
priori a formidable task, but progress can be made by recognizing that this is a semi- 



^As usual, X denotes the fraction of the hadron longitudinal momentum carried by a typical parton 
involved in the scattering: x ^ Q'^ /s, with a typical transferred momentum. Thus, the high energy 
limit (s — > 00 at fixed Q^) is the same as the small-x limit. 
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classical regime, for which a classical effective theory can be written [ This is a theory 
for the hadron wavefunction in the infinite momentum frame, and it looks remarkably 
simple: It is just classical Yang-Mills theory in the presence of a random distribution of 
colour charges which propagate nearly at the speed of light, and whose internal dynamics 
is "frozen" by Lorentz time dilation. The classical gauge fields represent the "small-x 
gluons", i.e., the gluons with small longitudinal momenta, which are the relevant degrees 
of freedom for high-energy scattering. The fast moving charges represent the partons 
with relatively high longitudinal momenta, like the valence quarks, which are not directly 
involved in the scattering, but act as sources for the small-x gluons. At saturation, the 
colour fields are strong, and the classical theory must be solved exactly. 

Originally, this theory has been formulated as a model to study saturation effects in a 
large nucleus (A 1) at not so high energies [|^, |, ^, ||, |^. In that case, the gluon density 
is high because of the presence of many colour sources (the Ax valence quarks) already 
at "tree- level" . Subsequently, it has been shown that this effective theory is consistent 
with the QCD evolution towards higher energies [ TU, TT|: The evolution modifies the 



classical colour source by incorporating the quantum gluons with longitudinal momenta 
above the small-x scale of interest. This involves a renormalization group equation which 
can be seen as a functional non- linear generalization of the BFKL equation (see Sect. 3 
below) . 

We have dubbed the gluonic matter described by this effective theory a Colour Glass 
Condensate (CGC) [|ll|, HI: 



• ''Colour" since gluons are "coloured" under SU(3). 

• ''Glass " because of the analogy with the mathematical description of systems with 
frozen disorder, called "glasses". Here, the "frozen disorder" refers to the random 
distribution of time-independent colour charges, which is averaged over in the cal- 
culation of physical observables. 

• "Condensate " since, at saturation, the gluon modes have large occupation numbers, 
of order l/tts (corresponding to strong classical fields A ~ !/(?), which is the maxi- 
mal value allowed by the repulsive gluon interactions. This is a Bose condensate. 

It is my purpose in this talk to review the physical and mathematical foundations of the 
CGC picture, and describe some of its consequences for gluon saturation [ |13l H, geo- 



metric scaling in deep inelastic scattering [ [T^, [1^, , and the unitarization of hadronic 



cross-sections [ |TB[. I will also discuss "limiting fragmentation" as an incentive towards 
the renormalization group description of multiparticle production in heavy ion collisions. 
Other consequences and recent applications to the phenomenology of deep inelastic scat- 
tering and heavy ion collisions will be briefly mentioned towards the end. More detailed 



discussions can be found in specialized contributions to these proceedings [ 19 



2. FROM BFKL TO GLUON SATURATION 

Consider deep inelastic scattering, for simplicity. We are interested in the high-energy, 
or small-x, limit, with x ~ Q'^/s ^1. It is convenient to regard this process in a special 
frame (the dipole frame) in which most of the total energy is carried by the hadron, while 
the virtual photon 7* has just enough energy to dissociate before scattering into a quark- 
antiquark pair (a colour dipole), which then scatters off the gluon fields in the hadron. 
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(See Fig. |T|.) In this frame, the DIS cross-section can be related to the dipole-hadron 
scattering amphtude A/'r(r_|_) [ which to lowest order (two gluon exchange) is in turn 
proportional to the gluon distribution in the hadron wavefunction: 

Here, r^ = x^ — is the transverse size of the dipole (with the quark at x± and the 
antiquark at y±), r = ln(l/x) is the rapidity gap between 7* and the hadron (which is 
almost the same as the hadron rapidity in the considered frame) , and xG{x, Q"^) / nR^ is the 
number of gluons with longitudinal momentum fraction x and transverse size r± = 1/Q 
per unit rapidity per unit transverse area. Throughout, I assume that the dipole is 
sufficiently small to be perturbative: Q"^ = l/rj_':^ ^qcd- 




Figure 1. Deep inelastic scattering in the dipole frame 

In writing eq. (|2.1j ), I assumed the hadron to be homogeneous in the transverse plane; 
that is, the gluon density was taken to be the same at all the impact parameters within 
the hadron disk of radius R. (A more realistic impact parameter dependence will be 
considered in Sect. 7.) Note that the scattering amplitude (^]l|) is not sensitive to the 
details of the longitudinal distribution of gluons inside the hadron, but only to their 
projected distribution in the transverse plane. Physically, it is so since the low-energy 
virtual photon has a large longitudinal wavelength oc so it scatters coherently off all 
the partons in a longitudinal tube of transverse size rj_. 

The F2 data at HERA show a rapid growth of the gluon distribution with 1/x, in 
qualitative agreement with the predictions of perturbative QCD to which I now turn. 

2.1. The small— X problem of the BFKL approximation 

Within perturbative QCD, the enhancement of the gluon distribution at small x pro- 
ceeds via the gluon cascades depicted in Fig. |^. Fig. |[a shows the direct emission of a 
soft gluon with longitudinal momentum^ = xP^ <^ by a fast moving parton (say, 

use light-cone vector notations: = (A:° ± k^)/\/2, k± = {k^,k^). The total 4-monientum of the 
hadron reads = with large P'^. 
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a valence quark) with = xoP~^ and 1 > xq x. Fig. Q.b displays the lowest-order ra- 
diative correction which is of the order (with A^^ the number of colours, and ctg = C(sNc/t^) 




a) 




Figure 2. a) Small-x gluon emission by a fast parton; b) the lowest-order radiative 
correction; c) a gluon cascade. 
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relative to the tree- level process in Fig. ^a. This correction is enhanced by the large 
rapidity interval At = ln(a;o/a;) available for the emission of the additional gluon. A 
similar enhancement holds for the gluon cascade in Fig. |^.c, in which the succesive gluons 
are strongly ordered in longitudinal momenta: p~^ ^ pf P2 ^ ■ ■ ■ ^ ^ A;"'". 
This gives a contribution of relative order ^ (^a^ In . Clearly, when x is so small 
that ln(a;o/a;) ~ l/^s, all such quantum "corrections" become of order one, and must be 
resummed for consistency. This gives the number of gluons with longitudinal momentum 
fraction x (or rapidity r = ln(l/x)) within this "leading-log" approximation: 

^ = xG{x,Q^) ~ e""^" = (2.3) 
dr 

with uj a pure number. We have tacitly assumed that all the gluons in the cascade 
have transverse momenta of the same order, namely of order Q. A more refined treat- 
ment, based on the BFKL equation [ 0], allows one to compute uj and specifies the 
Q^-dependence of the distribution. One obtains: 



xG{x, 



where u 



BFKL 



A^Q"^ exp < LuasT 



1 



In 



A2 



(2.4) 



41n2 2.77, /3 = 28^(3) ~ 33.67, and A is a reference scale, of order Aqcd- 
Eq. ( p^.4D exhibits two essential features of the gluon distribution in the BFKL approxi- 
mation: (a) its exponential increases with r, and (fe) its diffusive behaviour in InlQ"^ / A'^) , 
with diffusion "time" equal to r. These features are responsible for the main difficulties 
of this approximation in the high energy limit : 

(a) Violation of the unitarity hound : At high energy, hadronic cross-sections are pro- 
portional to the gluon distribution (see, e.g., eq. (pTTP). Then, eq. ( p^) predicts 
crtot(s) ~ s"^"", which violates the Froissart bound 
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(b) Infrared diffusion : With increasing s, the typical transverse momenta carried by the 
gluons within the BFKL ladder diffuse into the non-perturbative region at < A2. 
This contradicts the use of perturbation theory. 

2.2. The saturation momentum 

When resumming the ladder diagrams as above, we have assumed that the emitted 
gluons propagate freely, that is, we have neglected their interactions with gluons at the 
same rapidity from other cascades (cf. Fig. |^). Such interactions are not enhanced by a 
large logarithm, but may nevertheless become important when the density of the available 
gluons is large enough. Specifically, the interaction probability for gluons from different 




Figure 3. Gluon recombination in the hadron wavefunction. 
parton cascades can be estimated as: 

.(Q') X n(.. 0^) ~ ^ X l^^> (2.5) 

where cr(Q^) is a typical cross-section for gluons with transverse size 1/Q, and n{x, Q"^) is 
the density of the gluons of a given colour, and with given x and l/Q, in the transverse 
plane. This probability becomes of order one for lower than a critical value : 

r)2/ \ ^ (^sNc xG{x,Ql) 

that we shall refer to as the saturation momentum. More precisely, this is the scale at 
which non-linear effects start to be important in the hadron wavefunction: 

For ^ QKt), the non-linear effects are negligible, and linear evolution equations 
(like BFKL or DGLAP) apply. In particular, one can estimate the saturation scale by 
inserting the BFKL approximation ( |2.4| ) into eq. (|2.6|) . This gives [ [1^ : 



Qlir) = AV"^" , c = [-(3 + ^(3{(3 + 8cu) ]/2 = 4.84... (2.7) 

For Q"^ ^ Qlij)i the non-linear effects are essential, and are expected to tame the 
growth of the gluon distribution with r. It has been predicted almost twenty years ago 
[ m that a "saturation" regime should be reached in which recombination equilibrates 
radiation (see also Refs. [ 0, |, |^). To verify this conjecture, a formalism is needed to 
perform calculations in the non-linear regime. This will be described in Sect. 3. 
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For a nucleus, xGa{x, Q'^) oc A and vri?^ oc A^^^, so eq. (PD predicts qc A^/^. To 
summarize (with (5 ~ 1 /3 and A ~ ca^ in a first approximation) : 

Ql{x,A) ~ (2.8) 

which shows that an efficient way to create a high-density environment is to combine large 
nucleai with moderately small values of x, as currently done at RHIC. Schematically: 

— ep at HERA ~ heavy ions at RHIC 

— eA at eRHlC ~ heavy ions at LHC 

Eq. (|2.8|) also shows that the saturation momentum becomes a hard scale {Ql{x,A) 3> 
A.qcd) sufficiently large nuclei and/or large enough energies. Once this condition is 
satisfied, we can reliably use weak coupling techniques to study the evolution towards 
even higher energies. This strongly suggests that the high energy limit of QCD should be 
within the reach of perturbation theory. 

2.3. Saturation vs. Unitarization 

Since unitarity violations by the BFKL approximation are related to the unlimited 
growth of the gluon distribution, we expect that, converserly, saturation effects should 
restore unitarity for hadronic cross-sections at high energies. Consider, e.g., the dipole- 
hadron scattering amplitude, for which the unitarity bound reads N'r{rj_) < 1. (The 
upper limit A/'r(r_|_) = 1 corresponds to "blacknesss" : the dipole is completely absorbed.) 
Of course, BFKL violates this bound for high enough energies. But a brief inspection of 
eqs. (|2.1|) and ( p.6|) reveals that the unitarity violations show up precisely at the scale at 
which BFKL is expected to break down because of non-linear effects: 

Ar,(rx) ~ 1 for r^r^l/Qsir). (2.9) 

This suggests that unitarization effects and saturation are different aspects of the same 
non-linear physics, and cannot be dissociated from each other. We shall see below that, 
indeed, the non-linear effects responsible for saturation cure the unitarity problem as well. 



3. THE EFFECTIVE THEORY FOR THE CGC 

According to eq. ( |2.6|) , the gluon distribution is large at saturation: xG(x, Q^) ~ 
for ^ (''")) which suggests the use of semi-classical methods. One can write a 
classical effective theory whose general structure is fixed by the kinematics of the infinite 
momentum frame [ ^ : The "fast partons" move nearly at the speed of light in the positive 
z (or positive x~^) direction, and generate a colour current = S^^pa- By Lorentz 
contraction, the support of the charge density pa is concentrated near z = t, or x^ = 0. 
By Lorentz time dilation, pa is independent of the light-cone time x~^. 

The classical equation of motion (EOM) reads therefore: 

(D^F-'naix) = 5^+Pa(x-,Xx). (3.1) 



Physical observables are obtained by averaging the solution to this equation over all the 
configurations of p, with a gauge-invariant weight function Wr[p] which depends upon 
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the dynamics of the fast modes. For instance, the dipole-hadron scattering amphtude is 
computed as (in the eikonal approximation [ ^ ) : 



tT{v\x^)V{y^)))^ = jDpWr[p]tT{v:^^[p]VyAp]), (3.2) 

where V'^{x±) and V{yj_) are Wilson hnes describing the interactions between the fast 
moving quark, or antiquark, from the dipole and the colour field in the hadron: 

0(xx) = Pexpji^ J dx- A+{x-,x±)T''Y (3.3) 

and ^^[p] is the solution to the EOM ( p.l| ) in the covariant gauge: —V\A'^ = p. The 
average over p in eq. ( |3.2| ) is similar to that performed for a spin glass [ [12|. 

As already mentioned, this "glassy" description is an effective theory, that is, it holds for 
a given (small) value of x, i.e., for a given longitudinal momentum = xP^ <^ . The 
source pa is the colour charge density of the "fast partons" (the valence quarks and the 
gluons with momenta p"^ ^ in the parton cascades in Fig. |^). The classical solution 
A[p] represents the last gluon (with momentum fc"*") in these cascades. The non-linear 
effects in the classical EOM describe gluon recombination at the soft scale (see Fig. 

The corresponding effects at momenta p"^ ^ are included in the definition of p, 
that is, in the weight function WV[p]- 

All the dependence upon the separation scale A;"*" = xP^ is carried by the weight func- 
tion (via r = ln(l/x)) : When x is further decreased, say to xi < x, new quantum modes 
are effectively "frozen" (those with longitudinal momenta XiP^ < < xP'^), and must 
be included in the effective source at the lower scale XiP~^. This can be done via a one- 
loop background field calculation, and leads to a renormalization group equation for Wr[p] 
which shows how the correlations of p change with increasing r [ |10|, |ll|. Schematically: 



which describes diffusion in the functional space spanned by pa{x~ , x±) [ |2^. The kernel 
x[p], which plays the role of the "diffusion coefficient", is positive definite and non-linear 
in p to all orders (see Refs. [ |TT|, |T2[ for an explicit expression and more details). 

Eq. ( p.4| ) resums both the large energy logarithms, i.e., the terms ~ (a^lnl/a;)", and 
the leading high density effects, i.e., the non-linear effects which become of order one at 



saturation. Approximate solutions to this equation have been obtained in Ref. [ |T^, 
and will be described in Sect. 5 below. An exact, but formal, solution in the form of 
a path-integral has been constructed in Ref. [ where the stochastic nature of this 
equation has been fully elucidated, and the equivalent Langevin equation (describing a 
random walk on a group manifold) has been written down. Both the path-integral and 
the Langevin formulation are well suited for numerical simulations on a lattice [ |25[| . 



Eq. (|3.4|) is a functional equation, but can be transformed into ordinary evolution 
equations for the physical observables of interest. For instance, by taking a derivative 
w.r.t. r in eq. (|3.2|) , and using eq. ( ^^ for dWr/dr, one can obtain an equation for the 
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Figure 4. Particle production in heavy ion collisions. 



evolution of the scattering amplitude with r [ |TT|. Note, however, that in general this 
is not a closed equation (the 2-point function of the Wilson lines is coupled to a 4-point 
function), but only the first step in an infinite hierarchy of coupled equations. A closed 
equation for N'r{x±, y±) can be nevertheless obtained in the large- A'^c limit. This is known 
as the Balitsky-Kovchegov (BK) equation [ ^ ^], and reads schematically]^: 



^ = OsKbekl ® a/; - asAfr ® Afr , (3.5) 
or 

that is, in addition to the linear BFKL terms, it contains also a quadratic term which 
enforces unitarization (see Sect. 7). 

So far, the picture of the quantum evolution has been developed fully in the hadron 
infinite momentum frame: In the language of Fig. [^, the dipole rapidity has been kept 
fixed, while the hadron has been accelerated to higher and higher energies, with the result 
that the hadronic wavefunction (described here as a CGC) evolves according to eq. (|3l^) . 
Alternatively, via a change the frame, one could use the increase in the total energy to 
accelerate the dipole, and study the evolution of its wavefunction with r. In the language of 
eq. (|3.2|), this would correspond to keeping unchanged the weight function WVo[p] (^o is the 
rapidity of the hadron which is now fixed), but modifying the scattering operator to allow 
for additional gluons in the (evolved) dipole wavefunction. Since physics is boost invariant, 
both descriptions should lead to the same evolution equations for physical observables like 
A/'r(x^, y±). And, indeed, the evolution equations for Wilson line correlators derived from 
the CGC [ |1T| (in particular, eq. (p.5|)) are identical to the equations originally obtained 
by Balitsky [ ^ and Kovchegov [ ^ from studies of the dipole evolution. It has been 
first recognized by Weigert [ EBI that Balitsky's hierarchy of coupled equations can be 



compactly summarized into a single functional equation, which is essentially equivalent 
to eq. (|3.4}). More recently, (See Ref. [ ^ for a comparison between the two approaches.) 



4. LIMITING FRAGMENTATION 

The renormalization group description presented above is based on the separation of 
scales in rapidity (In 1/x) for the degrees of freedom in the hadron wavefunction. This way 
to formulate the problem is certainly useful for DIS, where one can tune the kinematics 

similar non-linear equation was originally suggested in Ref. [ || and proved in [ |^] in the double-log 
approximation. More recently, Braun has reobtained this equation by resumming "fan" diagrams [ ^ . 
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Figure 5. Limiting fragmentation in rapidity distributions at RHIC. 



to observe partons with a given value of x. Here, I would like to give you an example 
which shows that the same strategy can be also useful for hadron-hadron collisions, where 
partons with all values of x are a priori involved in the scattering. 

Consider multiparticle production in, say, heavy ion collisions at RHIC. A typical ra- 
pidity distribution of the particles emerging from such a collision is shown in Fig. 4. The 
leading particles are shown by the solid line and are clustered around the projectile and 
target rapidities. The dashed line is the distribution of the produced particles. Remark- 
ably, this latter appears to reflect the rapidity distribution of the partons in the colliding 
hadron wavef unctions. In particular, it preserves the separation of scales in rapidity. 

The most striking evidence in this sense comes from the feature of the data known as 
"limiting fragmentation" : When plotted as functions of y—yproj, the rapidity distributions 
measured at RHIC at different energies are to a good approximation independent of 
energy (see Fig. 5; the final version of these data, including the run at 19.6 GeV, is by 
now available in Ref. [ ^). The meaning of this becomes clearer if one notices that 
Vproj ~y ^ lii(l/a;), with X the longitudinal momentum fraction of the produced particle. 
Then, Fig. 5 shows that, with increasing energy, the "fast" (large x) degrees of freedom 
do not change much, while new degrees of freedom appear at smaller values of x. In other 
terms, the rapidity distributions of the produced particles are functions of x alone, and not 
of the total energy: This is similar to the Bjorken scaling of the parton distributions, and 
is naturally explained by assuming that hadron interactions are short-ranged in rapidity, 
so that the hadron distributions in rapidity reproduce the corresponding distributions of 
the liberated partons. 
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5. GLUON SATURATION 



In order to solve the RGE (|3.4]) , one needs an initial condition at low energy (say, at 
Xq ~ 10~^ ■ • ■ 10^^). For a large nucleus, at least, it is reasonable to assume that the only 
colour charges at such a large x are the valence quarks, which moreover are not correlated 
with each other, because they typically belong to different nucleons [ § . Then, the initial 
weight function Wr^ [p] is simply a gaussian, with 2-point function [ ||, [1^ 



{pa{x±)pb{y±))o = SabS'^^\x_L - y±) Ho, po=2asA/R\, (5.1) 

which measures the colour charge density of the valence quarks in the transverse plane. 

When increasing r, the quantum evolution will modify this simple initial condition by 
introducing correlations and non-linearities (that is, in general the evolved weight function 
will not be a gaussian any longer [ ^ H^)- The details of the evolution depend upon the 
structure of the kernel x sq. ( |3.4| ). For the present purposes, it suffices to say that x 
depends upon p via the Wilson line (|3.3| ). This allows for an intuitive distinction between 
the high-density and low-density regimes [ Q : 

— At low energy, or large transverse momenta kj_ 3> Q1{t), we are in a dilute regime 
where fields and sources are weak, and the Wilson lines can be expanded to lowest order: 
V'^{x±) ^ 1 + ig J dx~ A'^ {x^ , X ±) . Then, eq. (|3.4| ) reduces to the BFKL equation for the 
2-point function {p{k±)p{—k±))r = pT{k±) [[l^], with solution (cf. eq. (|2.4| )) : 



p,{k^) ^ ^/iofcie-"^^ (5.2) 

As compared to eq. ( |5.1| ), we note the emergence of transverse correlations, as well the 
rapid, exponential, increase with r. 

— At high energies, or low momenta k'^ ^ Qs(t), the colour fields are strong, ^ 1/g, 
so the Wilson lines rapidly oscillate and average away to zero: V ^V^ ^ 0. Then the 
kernel x becomes independent of p, and the r.h.s. of eq. (|3^) simplifies drastically [ |14| : 



Or 2tx Jk_L Spa{k_i)6pa{-k_i 
This is the standard diffusion equation for the Brownian motion. The corresponding 



2-point function increases only linearly with the evolution "time" r [ |T2| : 
Pr{k^) = {r- r{k^)) ^ ^ ^ In ^ , (5.4) 



TT -KCOis kj_ 

that is, logarithmically with the energy, since r ~ Ins. (In writing eq. (|5.4|), I have also 



used eq. (|2.7|) for the saturation scale). For most purposes, a logarithm is as good as a 
constant. We conclude that, at low momenta k±_ -C Qs{t), the colour sources saturate, 
because of the non-linear effects in the quantum evolution. 

A second important feature of eq. ( p. 4] ) is that it vanishes like k'j_ when k"^ ^ 0. This 
means that the saturated gluons screen each other in such a way that colour neutrality 
is achieved over a typical transverse area 1/Q1{t). This has important consequences 
for the infrared behaviour of the perturbation theory, and also for the unitarization of 
cross-sections at high energy [ (see Sect. 7 below). 
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Figure 6. (a) Gluon phase-space density as a function of k±; note the change in behaviour 
below Qg. (b) The same as (a), but for two values of x. 



Gluon saturation can be seen also in the gluon distribution at low k± [ |T3|, |T4| . Let 



1 d'N _ 1 dxG{x,kl) 
~ TiR^ drd^k^ ~ TiR^ d^k^ ^^-^^ 

denote the gluon density in the transverse phase-space. We have (t)r{kj_) oc fiT{k±)/k'j_ 
[ 0]. This relation, together with the previous results for the charge-charge correlator 
HT{k±), implies the graphical representations in Fig. 6, which are interpreted as follows: 
At very large k± ^ Qs{t), one finds the expected bremsstrahlung spectrum ~ l/A;^, 
due to radiation from independent colour sources (cf. eq. (|5Tl|)). For lower momenta, but 
still well above Qs{t), the spectrum gets softer ~ thanks to the BFKL evolution 

(cf. eq. ( |5.2| )), but the density increases rapidly with the energy, as x~^ with A = uag- 
Finally, at k± ^ there is marginal saturation (in the sense of only a logarithmic 

increase) with respect to both l/k'j_ and 1/x (cf. eq. ( ^.4|) ). Thus, with decreasing x, the 
new gluons are predominantly produced at large transverse momenta ^ Qs{t), where the 
density is lower, and the repulsive interactions are less important. 

6. GEOMETRIC SCALING 

At saturation, the gluon phase-space density defined in eq. ( ^.5|) reads (see Fig. 6) : 



In addition to its saturation features alluded to before, and to the l/ag enhancement, 
this function exhibits another interesting property, generally referred to as "geometric 
scaling": It depends upon the two kinematical variables r and k'j_ only via the ratio 
k']_/Ql{T) (the "scaling variable"). This is a natural property at saturation, where there 
is only one intrinsic scale, the saturation momentum. Other quantities, like the dipole- 
hadron scattering amplitude A/'r(r^), show in this regime a similar behaviour [ |T4| . 
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Figure 7. HERA data on the cross section for 7*p DIS from the region x < 0.01 and 
Q'^ < 400GeV^ plotted versus the scahng variable T = Q'^/Ql{x). From Ref. [ |15|. 



But at higher momenta kj_ 3> QI{t), or shorter distances r± <C 1/Qs{t), where satura- 
tion effects should not be important, there is a priori no reason to expect such a scaling. 
And indeed, the gluon distribution in eqs. ( |2.4| ) or ( [5.2| ), or the scattering amplitude in 
the BFKL approximation (which follows from eqs. ( p.l|) and (|2.4| )) : 



1- 91.^(,„91V 



(with the notation = 1/r^) show no obvious scaling. 

It may thus appear as a surprise that the HERA data on the total cross section for 7*p 
deep inelastic scattering show nevertheless geometric scaling to a quite good accuracy for 
X < 0.01 and all up to relatively high values ~ 400GeV^ [ (see Fig. ^. Such 



Q are significantly higher than the estimated value of the saturation scale at HERA, as 
extracted from fits to F2 within the "saturation model" [ ^ : — 1 • • • 2 GeV^. 



Motivated by this experimental observation, we have reconsidered the scaling properties 
of the dipole-hadron scattering amplitude at momenta above the saturation scale [ 
[TTH , and discovered that the BFKL solution shows approximate scaling within a window 
< ~ <5f/A^. The upper limit is of O(lOOGeV^), in agreement with the data. 

The basic argument is simple enough to be explained here: Let me replace A^ by Q1{t) 
as the reference scale in eq. (|6.2|), with QI{t) given by eq. (|2.7|): 

ln(QVA2) = ln(QVg^(^)) ^ ^-^^ ^g g^ 

This immediately yields, with 

Ql{r) 




AA,(g^) ~ exp|-A,ln-^-— ^(In-^l \. (6.4) 
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Note that the exponent in (|6^) vanishes for = Qlij)^ as expected (cf. eq. ( p.9|) ). 
From the equation above, we see that scahng emerges provided is close enough to 
Qlij) (ahhough still above it) for the second term in the exponent to be negligible. This 
provides the upper limit on alluded to above. 

To conclude this section, let me mention two recent applications of the idea of geometric 
scaling at saturation to the phenomenology of heavy ion collisions at RHIC. In Ref. [ 



it has been argued that the RHIC data for the transverse momentum distributions of the 
produced hadrons show geometric scaling in their dependences upon rriT and centrality: 
they are universal functions of itit/ ^s{h), with h = the impact parameter. In Refs. [ 
0, |3^, it has been argued that the scaling violation via the running of the coupling 
constant [a^ — > asi^Q^r)) in eq. ( |6.1| )] may explain the observed centrality dependence of 
the multiplicity at RHIC. 

7. FROISSART BOUND FOR DIPOLE-HADRON SCATTERING 

Let me now return to the dipole-hadron collision, and address the fundamental question 
of the asymptotic behaviour of the total cross-section as s — cxo. Since this is obtained 
by integrating the scattering amplitude J\fr{x±, y±) = Afr{r±, b±) (with r± = x±— y± and 
bi_ = {x± + y±)/2) over all the impact parameters b± (see Fig. 8) : 

atot(r,rx) = 2 J d%^Mr{r^,h^) , (7.1) 

it is clear that there are actually two issues that are concerned by this question: 
(i) the increase of the scattering amplitude with r at fixed impact parameter, and (m) the 
decrease of the scattering amplitude with 6_l at large impact parameters. The first issue 
is that of unitarity : the local scattering amplitude cannot exceed the unitarity bound 
■^T{f±,bj_) < 1. The second issue is that of confinement this is why hadron interactions 
have only a finite range. 

{i) The unitarization problem can be addressed within the perturbative framework. 
Clearly, the BFKL equation violates unitarity, but this is restored by saturation, which 
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Figure 9. The scattering amplitude as a function of rj_ for two values of r. 

ensures that Afr{r±,b±) ~ 1 for any r± ^ 1/<5s(t, Here, Qs{T,b±) is the local satu- 
ration scale, defined as in eq. (p^), but with xG{x, Q^)/nR'^ replaced by xG{x, Q^, b±) = 
dN /dT(Pb± (the gluon density at the given b±). 

The behaviour of Mr^r^, b±) as a function of r± at fixed b± is illustrated in Fig. 9 for 
two different energies. One sees a change in regime around r± ~ 1/Qs{t, b±) from "colour 
transparency" at small r± (cf. eq. ( p.l| )) to "blackness" at large r±. This behaviour can 



be obtained by numerically solving the BK equation (|3.5|) [ pO, 35, 36|. Some limiting 



regimes can be also studied analytically. For r± <^ 1/Qs{t, b±), one obtains [ 

AfAr^,b^) l-exp|-rl^^xG(a;,l/ri,6^)}, (7.2) 

which exhibits the unitarizing role of the multiple scattering (compare to eq. ( |2.1| ), which 
describes a single scattering). For r± ^ 1/Qs(t, &_l) one rather has [ 0, Q : 

K{r±M ^ 1 - exp|-l(lnQ2(^,5^)^y'|, (7.3) 

which shows geometric scaling with the scale set by the local saturation scale Ql{T,b±). 

With increasing r at fixed the saturation scale Q^(t, 6_l) increases and eventually 
becomes larger than the external resolution scale = l/r'j_ (see Fig. 8). Once this 
happens, the local scattering amplitude reaches the unitarity limit: N'r{r^,bj_) ~ 1 (i.e., 
the hadron becomes locally "black" ), and then remains constant when the energy is further 
increased. Since the gluon density is larger at the center of the hadron, the black disk 
first appears at b±_ = 0, and then extends radially with increasing r (see below). This 
transverse expansion explains why the total cross-section keeps growing even at very high 
energies. Let R{t, Q"^) denote the black disk radius. Then, by definition : 

K{Q^,b^)^l for b < R{t,Q^) (7.4) 

[with AfriQ^, bj_) = A/'r(r_L = 1/(5, b±) from now on], or, equivalently, 

Ql{r,b^) = for b = R{t,Q^). (7.5) 

In order to compute R{t,Q'^), one needs to understand the evolution of the scattering 
amplitude with r at impact parameters outside the black disk, i.e., in the "grey area" at 
b > R{t, Q'^) (see Fig. 8). This brings us to the problem of confinement. 
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(ii) If 6 > R{t,Q^), the gluon density is low and J\fr{Q^, bj_) is simply proportional 
to the local distribution function xG{x,Q'^,b±), cf. eq. ( |2.1| ). That is, the dipole is 
probing the gluon fields with momenta < created at bj_ by all the colour sources 
in the hadron. Because of confinement, such fields must fall off exponentially when the 
transverse separation between bj_ and the colour sources is larger than the pion Compton 
wavelength l/rriT, (since is the lowest mass gap in QCD). We thus expect [ ^ : 

K{Q\ &±) oc e-2'"-'' for b - R{t, Q^) > l/2m, , (7.6) 

where twice the pion mass enters the exponent because of isospin conservation: The long- 
range scattering is controlled by pion exchange, and pions have isospin one, while gluons 
have isospin zero, so the scattering proceeds via the exchange of (at least) two pions. A 
similar exponential decrease with holds for the saturation scale [ (see Fig. 8). 

Since MriQ"^, b±) is rapidly decreasing at b± 3> Rij^Q"^), the total cross-section is 
dominated by the black disk: 

atot(r,g') ~ 2rf(r,g2). (7.7) 

Thus, the relevant question is: How fast expands the black disk with r ? 

As shown in Ref. [ |18[) and I shall briefly explain here, one can answer this question by 
combining the perturbative (but non-linear) evolution with r with the non-perturbative 
boundary condition ( |7.6| ) due to confinement. Perturbation theory applies since the points 
at which the hadron turns from "grey" to "black" (see Fig. 8) are characterized by 
sufficiently large gluon densities, as we shall see. 

Specifically, we are in a weak coupling regime provided the local saturation scale is hard: 
Qlij-i b) ^ ^%CD- Given eq. ( |7.5| ) and the fact that the dipole is small (Q^ ^ Is?q(j^), it is 
clear that there exists a grey corona at -R(r, Q^) < b < Rh{t) within which perturbation 
theory applies (since Ag^^ -C Qlir.b) <^ at any 6j_ within this corona). Here, RniT) 
is the radial distance at which the saturation scale falls down to Aqcd (see Fig. 8). We 
shall check a posteriori that this perturbative corona is wide enough to allow for the 
calculation of the expansion rate of the black disk. 

For impact parameters within this corona, MriQ'^.b^) <^ 1, so one may expect the 
linear, BFKL, approximation to apply. This turns out to be correct eventually, but the 
complete argument is more subtle [ |18[: By itself, the BFKL equation would allow for 
long-range interactions between the incoming dipole and the colour sources in the hadron. 



but in the full equation ( p.5|) such interactions are suppressed by the non-linear effects 
whenever the transverse separation exceeds l/Qs{T,b±). Physically, this reflects the fact 
that the saturated gluons are colour neutral (cf. eq. ( p.4| )), and therefore couple to the 
incoming qq pair only via dipole-dipole forces, which have a rapid fall-off with bj_, and 
therefore do not contribute significantly to the scattering in the grey corona. Rather, 
the dominant contributions come from the nearby colour sources, i.e., the sources located 
within a saturation disk around the impact parameter b± (see Fig. 8): 

< Q;\T,b^). (7.8) 

To describe such short-range contributions, one can indeed rely on the BFKL equation, 
but which is supplemented with an infrared cutoff ~ Qs{t, b±), to remove the unphysical 
long-range contributions. 
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The crucial point is that the long-range interactions are cut off by saturation effects 
already at the short scale l/Qs{T,b±), which is much shorter than the soft scale 1/Aqcd 
at which the confinement plays a role. Because of that, the perturbative evolution with r 
is not sensitive to the transverse inhomogeneity in the hadron, which manifests itself only 
on the very large scale l/m^r ~ 1/Aqcd (cf- eq. ( [7.6|) ). That is, the evolution proceeds 
quasi-locally in the impact parameter space, in such a way that the 6^-dependence of the 
scattering amplitude factorizes out, and is fixed by the initial conditions. To summarize: 



grey ' ' BFKL 

~ exp <^ -2m^6 + wa^r " 2 A2 I ' '^^ 



where the fo^-dependence comes from eq. (|7.6| ), and in writing the second line I have 
used the BFKL solution ( |6.2|) (with the diffusion term omitted, as appropriate at high 
energies). This applies only for points in the grey area {h > R{t, Q'^)), but it can be used 
to estimate the black disk radius from the condition ([7.41). One obtains: 



which implies the following result for the total cross-section at high energy [ IT 



(^tot{s,Q^) ^ ^ (^) In^s as (7.11) 



2 V 



This saturates the Froissart bound, with a proportionality coefficient which is universal 
(i.e., the same for all hadrons), and which reflects the combined role of perturbative and 
non-perturbative physics in controlling the asymptotic behaviour at high energy. 

At this point, one should recall that the Froissart bound atot < (Joints [ ^ is a 
consequence of general principles (unitarity, crossing, and analiticity), but does not rely on 
detailed dynamical information. So, in reality, this bound may very well be not saturated. 
It so happens, however, that the measured cross- sect ions (e.g., for pp and pp scattering) 
show a slow, but monotonous, increase with s, which can be reasonably well fitted by 
a In^ s behaviour [ |3^ . This suggests that the Froissart bound is actually saturated 
in nature, and the mechanism described above provides a physical picture for such a 
saturation, with a definite prediction for the scale (Tq. 

From eq. (|7.9|), one can compute also the radius Rh{t) where the saturation scale 
decreases to Aqcd (beyond which perturbation theory fails to apply). According to 
eqs. ( |7.4|) -( [7^ , this is the same as the black disk radius for a large dipole with Q"^ ~ ^%cd'- 

R„ir) . (7.12, 

One can now compute the radial extent of the grey corona within which perturbation 
theory is applicable: 

RH{r) - R{t,Q') ^ ^Inf^. (7.13) 
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This is independent of r, and much larger than l/m.,^ (because of the large logarithm 
ln(Q^/A^)), which demonstrates the consistency of the previous calculation: During a 
"time" interval a^Ar ~ 1 (the typical rapidity increment at high energy), the black disk 
expands from R{t,Q^) to R{t,Q'^) + Lj/2m.„, cf. eq. ( |7.1CI| ), which is much smaller than 
Rh{t), and therefore still in the region controlled by perturbation theory. That is, the 
expansion of the black disk proceeds within the perturbative corona for intervals At which 
are large enough to allow for a controlled calculation of the rate of this expansion 



On the other hand, one cannot rely on perturbation theory alone to follow the expansion 
of the black disk for arbitrarily large "time" intervals At. For instance, it makes physically 
no sense to compute in perturbation theory the gluon distribution or the scattering am- 
plitude at very large impact parameters bj_ ^ Rh{t), where formally Qlij^h) <^ ^%cd- 
But if one insists in doing so, one finds that, because of the lack of confinement, the per- 
turbative evolution (say, according to the BK equation (|3.5D) generates long-range gluon 
fields which replace, at sufficiently large distances, the exponential fall-off of the initial 
distribution by just a power-law fall-off. If one pushes the perturbative expansion until 
the black disk enters this power-law tail, then its expansion rate speeds up, and eventually 
violates unitarity [ ^ . 



Clearly, this violation is an artifact of pushing perturbation theory beyond its limits of 
validity. To generate a significant power-law tail, the perturbative fields must propagate 
over distances many times the pion Compton wavelength. Such long-range fields are 
certainly unphysical: in the real world, they are removed by confinement. 

The true difficulty of perturbation theory is that it cannot be used to generate the 
exponential tail at large distances. But this is to be expected: such a tail can arise only 
from confinement. In any case, this difficulty has no incidence on the previous calculation 
of the rate of the expansion, since this requires only a limited evolution in r which is driven 
by short-range interactions. The prediction ( |7.11| ) of this calculation can be extended^ to 
arbitrarily large r, although the perturbative evolution becomes meaningless eventually. 

To conclude, while perturbation theory alone appears to be sufficient to describe uni- 
tarization at fixed impact parameter, one still needs some information about the finite 
range of the strong interactions in order to be able to compute total cross-sections. This 
is reminiscent of an old argument by Heisenberg [ |39[ which combines unitarity and short- 
rangeness (as modelled by a Yukawa potential) to deduce cross-sections which saturate 
the Froissart bound. Fifty years later, our progress in understanding high energy strong 
interactions allows us to confirm Heisenberg's intuition, and identify short-rangeness with 
confinement, and unitarization with saturation. 

8. APPLICATIONS TO PHENOMENOLOGY 

I conclude this review with a succint enumeration of recent applications of the concept 
of CGC to the phenomenology of deep inelastic scattering at HERA, and of relativistic 
heavy ion collisions at RHIC (with perspectives for LHC). More detailed discussions and 
more references can be found in [ . But let me start with a few words of caution: 

(a) The theoretical analysis underlying this concept relies heavily on the smallness of the 
coupling constant, which in turn requires a rather hard saturation scale: Q1{t) ^ ^%cd- 



Indeed, at any r, one can repeat this calculation and obtain an expansion rate consistent with eq. ( 7.10 ) 
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But for the currently available energies this condition is only marginally satisfied. Indeed, 
one estimates that _ 1 ■ ■ ■ 2 GeV^ at both RHIC and HERA [ 0. 

(6) The general formalism, and most of the theoretical calculations, have been devel- 
oped so far only to lowest order in a^. But for realistic comparisons with phenomenology, 
a NLO formalism (at least) is necessary. The NLO corrections to the BFKL kernel be- 
came finally available [ HO], but it turned out that resummation techniques are necessary 



to render these corrections meaningful [ 41]. Very recently, these methods have been 



used to estimate the saturation scale from the NLO BFKL equation [ ^ (via the same 
strategy as already used at LO; cf. Sect. 2.2). As a result, the logarithmic derivative 
As = d\n[Ql{T) / Is?']/ dr appears to decrease from the LO estimate ca^ ~ 1 (cf. eq. (p.6|) ) 
to As — 0.30, which is a reasonable value for the HERA phenomenology [ A fully 
non-linear NLO formalism is still to be developed (see however [ 0). 

(c) For applications to heavy ion collisions, the present formalism provides only the 
initial conditions, i.e., the wavefunctions of the incoming nuclei. Based on this, various 
methods have been developed to compute the melting of the CGC in the early stages 



of the collision, and the associated gluon production [ ||, ^ But very little 

is known from first principles about the subsequent evolution of the liberated partons. 
Their initial spectra and distributions can, and most certainly will, be modified by final 
state interactions, fiow, or thermalization, and may further change in the process of 
hadronization, but these various processes are not yet under theoretical control. It is 
therefore difficult, if not impossible, to unambiguously identify at this stage the effects of 
the initial conditions (say, of saturation) in the measured particle yields and correlations. 

This being said, it is remarkable that there exists already a significant amount of data, 
both at HERA and at RHIC, which hint towards saturation and the CGC, or at least are 
consistent with this physical picture. 

Let me start with DIS at small-x, for which the theoretical calculations are better 
under control, and quantitative predictions can be made. It has been shown in Refs. [ 



31] that the HERA data for both inclusive and diffractive proton structure functions can 



be well accounted for by a phenomenological model which incorporates saturation. The 
success of this simple model for diffraction is particularly significant, since the diffractive 
cross-section is more sensible than the inclusive one to dipole fluctuations with a large 
transverse size r^, and thus to saturation effects in the proton wavefunction. The same 
model has motivated the search for geometrical scaling in the HERA data [ |T3[, with the 
remarkable results discussed in Sect. 6. As shown in Ref. [ it is possible to use the 
data on the differential cross-section for vector meson production in DIS (say, 7*]? — pp) 
to obtain information about the impact parameter dependence of the dipole scattering 
amplitude. The phenomenological analysis in Ref. [ Q indicates a reasonable amount of 
blackness in the proton at central impact parameters already for x ~ 10~^. 



Consider heavy ion collisions now. In Refs. [|3^, 0, the particle production in Au — Au 
collisions at RHIC has been flrst analyzed from the perspective of the CGC. Specifically, 
the measured hadron multiplicities have been compared, in so far as their dependences 
upon centrality, rapidity, and total energy are concerned, to the respective predictions of 
the CGC picture. The agreement is quite good, in fact, remarkably good, given the many 
sources of uncertainty that I have mentioned before. 

A systematic approach, aiming at the numerical resolution of the classical Yang-Mills 
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equations which describe the scattering of two CGCs, has been developed and successively 
refined in Refs. [ In this approach too, one finds a reasonable agreement between 
the centrality dependence of the gluons liberated in the collision (as obtained from lattice 
simulations) and that of the charged hadrons measured at RHIC. 

For proton- nucleus (pA) collisions, and for DIS, there are also analytic calculations of 
the gluon and particle production, which include all multiple rescatterings off the strong 
colour field (the CGC) of the nucleus. The inclusive gluon production cross-section has 
been computed for pA [ ||, ^ and DIS [ [ 



, and a similar calculation for AA has 
Still for pA collisions, one has calculated the cross-sections for 
1^, photons and dileptons 



been attempted in [ ^ 



the production of jets [ photons and dileptons [ The charm production from 
the CGC in peripheral heavy-ion collisions has been investigated in [ See also the 
closely related work in Ref. [ Instantons in the saturation environment have been 
considered in Ref. [ KH] . At a conceptual level, there has been significant progress towards 
understanding the role of inelastic scattering for thermalization, and the way how kinetic 
theory could emerge from the classical field dynamics that applies at the earliest stages [ 
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